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Abstract 

The  cross  section  for  the  scattering  of  electrons  by  hydrogen  atoms 
at  zero  energy  has  been  calculated  using  a  variational  principle  of  the  Hulthen- 
Kohn  type.  For  singlet  scattering  we  have  used  a  three  parameter  trial'  function 
which  reduces  tc  an  H  wave  ftinction  at  small  electron  separations  and  which  has 
the  correct  asymptotic  form.  To  calculate  triplet  scattering  we  use  an  antisym- 
metrized  version  of  the  same  trial  function.  The  effective  range  for  singlet 
scattering  is  determined  from  the  scattering  length  already  calculated  and  the 
known  value  of  the  electron  attachm.ent  energy  for  H  .   For  triplet  scattering 
we  use  an  energy-dependent  trial  function  to  determine  the  effective  range. 
Using  the  calculated  values  of  the  above  parameters  we  find  the  cross  section 
as  a  fiinction  of  energy  for  the  range  0  -  0.02  volts.  The  values  of  the  cross 
section  are  then  used  to  calculate  the  resistivity  of  silicon  alloys  due  to 
neutral  impurity  Sv^^ttering. 
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1.  Introduction 

The  scattering  of  low-energy  electrons  by  atomic  systems  is  of 
some  interest  as  one  of  the  fundamental  processes  which  occur  in  ionospheric 
and  astrophysical  phenomena.  More  detailed  knowledge  of  the  scattering  pro- 
cess would  also  be  of  great  aid  in  xinderstanding  certain  aspects  of  gas 
discharges. 

The  simplest  case  of  this  sort  is  the  scattering  of  electrons  by 
hydrogen  atoms.  Measurements  of  the  cross  section  for  this  system  have  been 
made  by  BedersonL  -'  in  the  energy  range  above  2  volts.  However,  experiments 
for  the  energy  range  from  zero  electron  energy  to  thermal  energy,  a  range  of 
practical  interest,  are  extremely  difficult  to  carry  out.  At  present  there 
exists  no  data  for  this  energy  range. 

In  order  to  gain  some  knowledge  of  the  lovr-energy  scattering  it  is 
necessaiy  to  carry  out  a  theoretical  calculation.  The  problem  of  the  scatter- 
ing of  electtons  by  hydrogen  atoms  is  a  three-body  problem  and  thus  not  amenable 
to  exact  solution.  One  is  therefore  compelled  to  seek  approximate  methods  of 
solution. 

The  first  attempts  to  solve  the  problem  consisted  in  replacing  the 
three-body  problem  by  the  scattering  of  an  electron  by  a  center  force.  One 
thus  considered  the  scattering  of  electrons  by  a  Hartree  field  or  by  a  Fermi- 
Thomas  statistical  potential.  These  approximations  suffered  from  some  serious 
deficiencies.  For  one  thing  they  did  not  take  account  of  the  polarization  of 
the  atom  by  the  field  of  the  incident  electron,  nor  did  they  adequately  provide 
for  the  consideration  of  exchange  effects  which  are  extremely  important  for 
low-energy  scattering. 
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McDougall '-  -*  and  Chandrasekhar  and  Breen  ^  -•  have  treated  low- 
energy  electron-hydrogen  atom  scattering  in  the  one-body  approximation  using 
the  method  of  partial  waves  to  calculate  the  cross  section. 

The  need  for  more  accurate  approximation  methods  has  led  to  a  study 

of  the  applicability  of  variational  techniques  to  collision  problems. 

r22l 
Schwinger"-  -^  has  developed  very  powerful  techniques  which  are  based 

on  converting  the  Schrodinger  differential  equation  into  an  integral  equation 
involving  a  Green's  function  appropriate  to  the  particular  collision  problem 
under  consideration.  These  methods  give  stationary  expressions  for  the  scat- 
tered amplitude  arid  the  phase  shifts.  These  techniques  have  been  extended  by 
Horowitz  and  Friedman ^--^  to  handle  three-body  collisions. 

Other  variational  techniques  based  on  the  differential  equation  have 
teen  developed  by  Hulthen'-  -'.  Kohn"-  -•  has  developed  variational  techniques  for 
many-body  collision  problems.  Kato'-  -•  has  developed  a  variational  procedure 
which  gives  upper  and  lower  bovmds  on  the  phase  shifts  for  one-body  collisions. 
However,  methods  for  obtaining  bounds  for  the  phase  shifts  for  three-body  colli- 
sions are  still  lacking. 

The  variational  methods  have  been  applied  to  the  calculation  of  the 
electron-hydrogen  atom  scattering  cross  section  by  numerous  authors.  Vkssey  and 
Moiseiwitsch'-  -'*'-  -^  have  used  the  methods  of  Hulthen  and  Kohn  to  calculate 

ri2i 

the  elastic  scattering  cross  section  of  electrons  by  hydrogen  atoms.  Huang"-  -^ 
has  developed  his  own  technique  and  applied  it  to  electron-hydrogen  atom  colli- 
sions., taking  polarization  into  accoiint  but  neglecting  exchange  effects.  Massey 
and  Moiseiwitsch  have  shown  how  to  modify  Huang's  method  so  as  to  consider  ex- 
change effects.  Boyet  and  Horowitz'- J  have  applied  the  Schwinger  method'-  -'  to 
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calculate  the  direct  scattered  amplitude  for  elastic  collisions  of  electrons 
with  hydrogen  atoms.  Moiseiwitsch'-  -■  has  extended  Hulthsn's  procedure  to 
inelastic  processes  and  has  determined  the  phase  shdft  for  the  1S-2S  excita- 
tion of  hydrogen  atoms  by  electrons. 

None  of  the  above  calculations  are  well  suited  for  the  range  near 
zero  incident  electron  energy.  The  Hulthen-Kohn  variational  calculations  appear 
to  offer  a  better  approach  to  problems  in  this  energy  range  than  other  techniques. 
Born  approximation  procedures,  for  example,  are  good  only  for  high-energy  scat- 
tering, and  the  Schwini^er  variational  method  •- J  leads  to  integrals  which  cannot 
be  easily  evaluated  (cf.  [5j). 

In  a  similar  nuclear  scattering  problem,  namely  that  of  low-energy 

r22l 
neutron-proton  scattering.  Schwinger^  -•  has  shown  that  the  boimd  state  of  the 

deuteron  lying  near  the  continuiim  gives  the  dominant  contribution  to  the  cross 
section.  The  same  situation  should  hold  for  the  case  of  lovr-energy  electron- 
hjTirogen  atom  scattering,  since  there  is  a  bound  state  of  the  H  ion  lying 
about  three  quarters  of  a  volt  below  the  continuum^  -^  * '- -" .  By  using  a  varia- 
tional method  with  a  trial  fimction  suitable  for  the  H  ion  we  probably  in- 
clude all  of  the  polarization  effects  when  the  low-energy  electron  is  close 
to  the  iQ^drogen  atom.  The  bound  state  of  H  is  a  singlet  state,  so  that  we 
have  no  similar  guide  for  choosing  a  trial  function  for  the  triplet  scattering. 
However,  this  choice  is  not  crucial^ as  we  shall  show  that  the  singlet  scatter- 
ing is  the  dominant  featiire  in  the  process. 

It  is  the  pxirpose  of  this  paper  to  make  an  acciirate  calciilation  of 

the  cross  section  at  zero  incident  electron  energy.  We  shall  then  use  the  ef- 

[21  [3l 
fective  range  theory ^  j»l-'j  ^q   find  the  cross  section  at  small  but  finite 

electron  energies. 
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In  Section  2  we  define  the  phase  shifts  for  three-body  scattering 
by  using  the  S-niatrLx  formalism.  We  show  there  that  the  phase  shifts  are 
defined  only  if  the  space  part  of  the  wave  fiinction  is  symmetric  or  anti- 
symmetric in  the  coordinates  of  the  electrons.  We  then  extend  the  Kohn-Hulthen 
variational  principle  to  three-body  collisions.  Using  this  method  we  calculate 
the  singlet  and  tiriplet  scattering  lengths.  From  these  results  we  calculate 
the  total  cross  section  for  zero-energy  incident  electrons. 

In  Section  3  we  extend  the  effective  range  theory  to  three-body 
collisions.  Using  the  fact  that  a  bound  state  of  H"  exists  we  calculate  the 
effective  range  for  the  singlet  scattering.  We  next  introduce  an  energy- 
dependent  trial  function  and  use  it  to  find  the  effective  range  for  triplet 
scattering. 

With  these  results  we  calculate  the  total  cross  section  for  incident 
electron  energies  between  zero  and  thermal  energies. 

In  Section  h   we  compare  our  results  with  the  extrapolated  results  of 

ri9i 

Massey  and  Moiseiwitsch"-  -•  for  the  same  energy  range.  We  then  discuss  the 
sensitivity  of  the  cross  section  to  the  parameters  which  appear  in  the  effective 
range  theory.  Finally  we  obtain  an  estimate  of  the  accuracy  of  our  calculation. 

ErginsoyL  -■  has  used  the  scattering  of  electrons  by  hydrogen  atoms 
as  a  model  to  explain  the  effect  of  neutral  iit^jurities  on  the  resistivity  of 
semiconductors  at  low  temperatures.  In  Section  5  we  use  our  improved  values  of 
the  cross  section  to  calculate  the  resistivity  at  low  temperatures  and  compare 
these  with  the  values  obtained  by  Erginsoy  using  the  cross  sections  of  Massey 
and  Moiseiwitsch  •-  -• . 
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2.  The  variational  principle 

VJe  begin  by  considering  the  three-body  vrave  equation  for  S-vra.ve 
scattering  with  the  origin  of  coordinates  at  one  of  the  particles  whose  mass 
we  assume  infinite.  This  equation  is 

(1)  L  Kr^^g)  X  0 

where  L  is  the  operator 

here  r^  and  r_  are  the  distances  of  particles  1  and  2  from  the  origin,  nu  and 

m«  are  the  masses  of  the  electrons,  "fe  is  Planck's  constant  divided  by  2n,  V  is 

i?     2 

the  potential  energy  of  the  system,  ■«-—  V  is  the  kinetic  energy  operator  and 

E  is  the  total  energy  of  the  system.  For  the  case  of  scattering  of  electrons 
by  hydrogen  atons  (1)  becomes,  after  some  rearrangement  of  terms 

where  k  is  the  wave  number  of  the  incident  electron.  We  have  set  m,=  m«=  m  = 
the  mass  of  the  electron.  We  have  chosen  units  such  that  e  =  m  =  n  »  1,  where 
e  is  the  charge  of  the  electron.  In  these  \xnits  the  radivis  of  the  first  Bohr 
orbit  is  equal  to  \inity  and  the  binding   energy  of  an  electron  in  the  ground 
state  of  a  hydrogen  atom  is  E„  =  -  ■»  » 

The  meaning  of  the  other  quantities  may  be  seen  from  Figure  1. 

We  shall  now  rederive  the  Kato  principle  for  the  case  of  three-body 
scattering  and  apply  it  to  low-energy  electron-hydrogen  atom  scattering. 

Before  proceeding,  however,  it  will  be  necessary  to  generalize  the 
phase  shift  analysis  to  the  problem  of  three-body  scattering. 
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electron 


electron 


',z''2-'> 


proton  (infinite  mass) 


Figure  1 

For  two-body  scattering  with  central  forces  the  diagonal  elements  of 
the  S-roatrix  are  related  to  the  phase  shift  by 

For  this  case  the  spherical  harmonics  are  a  representation  in  which 
the  S-matrix  is  diagonal. 

For  the  case  of  two-body  scattering  with  tensor  forces  the  S-matrix 
is  no  longer  diagonal  in  this  representation.  In  order  to  define  the  phase 
shift  we  miist  therefore  diagonalize  the  S-matrix.  Schwinger  has  shown  for  the 
problem  of  neutron-proton  scattering  with  tensor  f6rces  that  the  eigenvalues 
of  the  S-matrix  are  related  to  the  phase  shifts  by  the  relation  given  in  ih) • 
We  can  define  a  phase  shift  generally  in  terms  of  the  diagonal  elements  of  the 
S-matrix  by  equation  (U).  In  the  case  of  the  scattering  of  electrons  by  hydrogen 
atoms,  the  spherical  harmonics  do  not  form  a  basis  for  diagonalizing  the  S-raatrix 
because  of  exchange  scattering. 
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In  order  to  diagonalize  the  S-niatrix  we  follow  the  method  due  to 
Kohn  who  has  shown  that  the  asymptotic  form  of  the  wave  equation  for  an  electron 
colliding  with  an  atomic  system  can  be  written  as 


(5) 


r(l^) 


^  ^±\\ 


,(P-)  e 


-ikr. 


ikr. 

-  &^  ^^ 

•^i    r. 


where  P.  is  the  angular  part  of  the  wave  function,  N.  is  a  normalization  con- 
stant, and  ^.   is  the  normalized  ground-state  wave  function  for  the  atom.  The 
index  i  mns  over  1  or  2  and  distinguishes  the  coordinates  of  each  electron. 
The  superscript  ^i   =  1,2,  refers  to  the  two  different  modes  of  collisions  which 
can  occxir,  depending  on  whether  electron  1  or  2  is  incident. 
The  S-matrix  is  now  defined  by  the  relation 

(6) 


^l  -  z:   3,,  a^ 


.  .  U.  .    • 


For  the  case  of  the  scattering  of  electrons  by  hydrogen  atoms  we  maj 
choose  two  f^ whose  asymptotic  forms  are 


(7) 


(8) 


r(l) 


Y 


r(l) 


(2) 


hh 


N2P2 


'-ikrn 


Ikr-i 


ikr-- 


-  (l+2ikf) 


,  0   (r^)  -  2ikg  - — 
Tn  /  ^o  2      ^   r. 


^o(--l' 


ikTj  ikr, 

«  (l+2ikf)  - 


2   .       ."'"2\  .^*^1 


jZl^(r^)  -  2ikg 


r.  ^0(^2) 


where  Y        is  the  wave  function  for  particle  1  incident  with  particle  2  bound 

(2) 
and  I        is  that  for  particle  2  incident  and  particle  1  bound.  0   (r^ )  and  J2J_(r-) 

are  the  normalized  ground  state  hydrogen  atom  wave  functions,  and  f  and  g  are 

the  direct  and  exchange  scattered  amplitudes. 
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The  elements  of  the  S-niatrix  are  clearly 


(10) 


1 


^11  "     ^22  "     -^  *  ^^^ 


^12   '     ^21  "     ^^^^       • 


we  must  now  find  the  linear  combinations  of  T^   and  Y^  '^  which  make 
the  S-matrix  diagonal. 

The  unitary  matrix  that  diagonalizes  the  S-matrix  is 

(u)  n     .  A.  (  ,    . 

The  particular  linear  combinations  which  diagonalize  the  S-matrix  are 
then  found  from  the  relation 

(12)  U  •  .  -     '  "  .1 


Thus  the  space  wave  functions  which  describe  the  scattering  must  be 
either  symmetric  or  antisymmetric  in  the  coordinates  r,  and  r^* 

The  eigenvalues  of  the  S-matrix  are  now  easily  found  by  using  the 
relation 
(13)  Sjj  =  US  u"^ 

and  are  given  by  relations 

2i6 
{Ih)  S^j^  -  e   -^  -  l+21k(f  +  g), 

2i62 
(15)  ^22  "  e  "  '    l+2ik(f  -  g)  . 
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From  the  above  discussion  vre  see  that  the  scattering  is  described 
by  two  phase  shifts  6,  ans  6_  corresponding  to  the  symmetric  and  antisymmetric 
wave  functions  respectively. 

V/e  now  proceed  to  the  generalization  of  the  Kato  variatiaaL principle. 
Consider  the  integral 

(16)        I 

where  L  =  E  -  H.  H  is  th3  Harailtonian  for  the  system  and  E  is  the  total  energy. 
Equation  (16)  can  now  be  written  as 


=  I  (uLv  -  vLu)  dr^  dtg  , 


(17)      I  =     [a(vj+  7^)v  -  v(vj+  V^)uJ  d-V^  dr^ 


|u7^v  -  v7^ujd  r^d  r^  +    j    Lv^v  -  vv^ujd  r^d  r^    , 


The  application  of  Green's  theorem  gives 


(18) 


I   fuV^v  -  vV^uldS^dr^  +  )   fuVgV  -  vVgiTldSgd  r^ 


where  V^  and  Vp  are  the  usual  gradient  operators.  dS-  and  dSp  are  the  surface 
elements  of  spheres  in  r, -ajid  r^-space  at  infinity. 

Let  us  now  choose  u  =  Y  to  be  a  solution  of  (1),  that  is,  LT  =  0,  and 
we  choose  v  =  T,  •  1.    is  chosen  so  as  to  have  the  same  asymptotic  form  as  ¥ 
except  that  it  may  differ  in  phase. 

VJ'e  take  the  asymptotic  forms  for  Y  and  f .  as 

(19)     Y  ->  —^  p-  Uos(kr3^+  e)  +  XQsin(kr^+  ®^o^^2^  -  ^1^  ^2  ' 


(20)      Y^  -> 


■^  1.   [cosdcr^-.  e)  +  Xg^sin(kr^.  e)]^o^^2^  1  ^1^  '"a  ' 
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where 

cot(6  -  e) 
(21) 


=  (cot  5  -  e)  i 

here  6  is  the  phase  shift.  From  the  previous  discussion  we  see  that  two  phase 
shifts  are  defined:  one  for  the  symmetric  wave  function  and  one  for  the  anti- 
symmetric wave  function,  so  that  we  shall  have  a  variational  principle  for 
each  phase  shift. 

If  we  now  insert  equations  (19)  and  (20)  into  (18),  and  use  (16), 
we  obtain  the  result 


If  we  now  let  w  =  Y  -  Y. ,  equation  (22)  becomes 


(22)         k(XQ-  Xq)  -  I  Y  L  Y^  d^'^drg  . 


-  I  Y^LY^  d  r^d  rg  +  j  wLw  d  r^c 


(23)         kXg  =  kXg  -    Y^LY^dr^drg*  I  wLw  d  r^d  r^ 

or 

(2h)     k  cot(5-e)  =  k  cot(6^-  e)  -  j  Y^LY^  d  r^d  r^  +  J  wLw  d ^^d  r^  . 

The  last  term  on  the  right  is  of  second  order  in  the  error  Y  -  Y.  ,  so 
that  (2li)  is  a  variational  principle  for  k  cot(6  -  ©) .  If  we  now  set  6=0, 
then  (2U)  reduces  to  the  Kohn-Hulthen  principle 

(25)         k  cot  5  =  k  cot  6^  -  I  Y^LY^  d  r^d  r^ 

which  is  the  form  of  the  variational  principle  we  shall  use  in  our  calculations. 

We  shall  now  apply  (25)  to  the  calculation  of  k  cot  6  for  the  scatter- 
ing of  zero-energy  electrons  by  hydrogen  atoms.  Equation  (3)  then  reduces  to 


(«'       pi*  '^  ^(^  *  4  -  4  -  ^  ^'^1-'^ 


2)  -  0  . 
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The  trial  functions  which  we  choose  for  use  in  the  variational  principle 
must  be  chosen  so  as  to  have  the  correct  asymptotic  form  as  given  by  equation  (20). 
Furthermore,  in  order  to  define  the  phase  shifts  the  trial  functionsimist  be  sym- 
metric or  antisymmetric  in  their  space  coordinates  r,  and  Tp. 

We  point  out  here  that  the  trial  functions  used  by  Massey  and 
Moiseiwitsch*-  -•  in  their  one-body  and  polarization  approximations  and  the  trial 
function  used  by  Huang  do  not  have  the  proper  symmetry^  so  that  the  phase  shifts 
can  not  be  defined*  The  use  of  these  trial  functions  for  the  calculation  of  the 
pliase  shifts  therefore  leads  to  meaningless  results. 

Since  we  are  dealing  with  electrons,  the  total  wave  function  (space 
and  spin)  must  be  antisymmetric  in  the  space  and  spin  coordinates  of  both 
electrons  in  order  to  satisfy  the  Pauli  principle.  For  the  case  of  the  sym- 
metric space  wave  function  the  electron  spins  must  be  aligned  antiparallel 
(singlet  scattering),  while  for  the  antisymmetric  space  wave  function  the 
electron  spins  will  be  aligned  parallel  (triplet  scattering). 

From  the  previous  discussion  we  can  tnerefore  define  two  phase  shifts, 
one  for  the  singlet  and  one  for  the  triplet  scattering. 

It  is  well  known  that  a  bound  state  of  H"  exists.  In  the  H  ion  the 
electrons  have  their  spins  aligned  antiparallel  (singlet  state),  Chandrasekhar'--' 
has  used  the  Ritz  method  with  a  trial  function  of  the  form 

e    e    +  e    e 

to  calculate  the  binding  energy  of  the  second  electron  in  the  ionj  his  result  is 
0.7019  volts. 

In  choosing  a  trial  function  for  our  phase  shift  calculation  for  the 
singlet  scattering,  we  should  like  it  to  reduce  to  the  trial  function  of 


^t  = 

Ae       -^e       "^  +  6       "^e       '^ 
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Chandrasekhar »- J  when  both  the  electrons  are  close  to  the  nucleus.  When  the 
incident  electron  is  far  from  the  hydrogen  atom  we  should  like  the  trial  func- 
tion to  have  the  correct  asymptotic  form.  For  the  case  under  consideration 
this  will  be  the  product  of  a  normalized  ground  state  hydrogen  atom  wave  func- 
tion with  the  sum  of  an  incident  plane  wave  and  a  scattered  wave.  With  these 
considerations  in  mind  we  choose  the  following  trial  function: 

(27) 

1—  VT-  21 

where  the  factor  containing  A  is  called  ^  and  the  remainder  is  1l,     A,  B  and  y 
are  the  parameters  tc  be  varied  and  1/  ^5"  n  is  a  normalization  constant.  In  the 
present  calculation  we  shall  assume  a  and  p  are  constant  and  use  the  values  ob- 
tained by  Chandrasekhar,  namely,  a  =  1.03925  and  p  =  0.28309. 

We  note  that  in  order  to  get  the  most  acc\irate  results,  a  and  p  shoiild 
also  be  varied.  However,  if  this  is  done  the  numerical  calculations  involved 
become  too  burdensome. 

The  symmetric  form  of  (2?)  will  describe  the  singlet  scattering.  For 
the  triplet  scattering  we  shall  use  the  antisymmetric  form  of  (2?). 

If  we  substitute  equations  (26)  and  (27)  into  (25),  it  takes  the  form 

(28)  X  «   lim   (k  cot  6)  =  -  a  a2+  b  B^+  c  AB  +  d  +  e  A  +  (f  -l)B 

o         y.  ^Q  Looo     oo^oj 

where  B  =  k  cot  5. . 

We  must  now  make  (28)  an  extremum  with  respect  to  A,  B  and  y» 
The  conditions  3X  /3A  =  aX  /dB  =  0  yield  the  relations 

(29)  2aA  +  cB  +  e     =     0 

000 

(30)  cA  +  2bB  +   (f-l)      =     0     , 

000  ' 
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which  give  for  A  and  B 

c  (f  -  1)  -  2b  e 

(31)  A     =      -^-2 0-0 

Ua  b  -  c"^ 
0   0        o 


c  e  -  2a  (f  -  1) 


c•5o^                      T,     -        0  0        0     0 
(.32;  B     =      2 

ha  b  -  c 
o  o      o 

The  quantities  a  ,  b  ,  etc.  are  functions  of  y>  which  are  obtained 
by  evaluating  the  integrals  which  occur  when  we  substitute  equations  (26)  and 

(27)  into  equation  (25).  This  involves  the  evaluation  of  a  great  many  integrals, 
most  of  which  are  elementary  in  nature.  These  calculations,  straightforward  but 
quite  lengthy,  are  given  in  Appendix  I.  Once  the  quantities  a  ,  b  ,...  are  ob- 
tained as  functions  of  y  w®  can  calculate  X  as  a  function  of  y  numerically,  sub- 
ject to  the  conditions  (31)  and  (32),  and  look  for  the  value  of  y  which  makes 

(28)  an  extremura. 

The  numerical  calculations  involved  were  carried  out  with  the  aid  of 
the  Burroughs  ElOl  digital  computer.  The  results  are  summarized  in  Table  Is 


Table 

I 

Singlet 

X 

0 

-0.1226 

Y 
0.12 

A 

o.oUih 

B 

-0.1282 

Triplet 

-0.^251 

0.1x6 

0.005^7 

-0.^256 

The  trial  wave  function  (27)  probably  does  not  take  adequate  account  of  the 
polarisation  of  the  hydrogen  atom  in  the  field  of  the  incident  electron.  In  order 
to  do  this  more  accurately  we  shall  modify  our  trial  function  by  making  it  de- 
pend explicitly  on  r^_,  the  interelectronic  distance.  Such  a  wave  function 
can  be  written  as 
(33)  T^  -  (1  +  Cj^r^g^  1  *  )^  » 


-  lU  - 


where  c,  =   0.3121h  and  is  the  value  found  by  Chandrasekhar.  I'he  values  of 
o  and  p  for  this  wave  function  are  a  =  1.07U78  and  p  =  O.I47758. 

The  quantities  c,  p,  and  c,  will  again  be  taken  as  constants  for 
this  calculation. 

The  integrals  involved  in  the  determination  of  a  ,  b  in  (28)  become 
somewhat  more  complicated  due  to  the  explicit  dependence  of  the  trial  function 
on  r-p.  Their  evaluation,  however,  is  still  fairly  straightforward  (see  Appendix  II)  ( 

It  shqjild  be  pointed  out  here  that  the  antisymmetric  form  of  equations 
(27)  and  (33)  probably  is  not  a  very  good  trial  function  for  the  triplet  scatter- 
ing. In  choosing  the  singlet  trial  function  we  were  guided  in  our  cl|pice  by 
the  fact  that  a  botind  state  of  K~  exists.  For  the  triplet  scattering  we  have 
no  such  guide.  There  is  therefore  no  reason  to  believe  that  antisymraetrizing 
the  singlet  trial  function  will  lead  to  a  good  trial  function  of  the  triplet 
scattering  when  the  electrons  are  close  together. 

We  shall  show  in  Section  k   that  the  singlet  scattering  gives  the  main 
contribution  to  the  cross  section,  so  that  eirrors  in  the  triplet  trial  wave 
function  will  make  very  little  difference. 

We  have  therefore  not  calciilated  X  for  the  triplet  scattering  using 
(33) •  F\irthermore ,  Massey  and  Moiseiwitsch  indicate  in  [I9j  that  the  polariza- 
tion correction  does  not  affect  the  triplet  scattering  to  any  great  extent. 

The  results  of  the  aforementioned  calculation  are  given  in  Table  II: 

Table  II 

X  Y  A  B 

o  ' 

Singlet  -0.1291  0.07  0.0308  -O.II472 
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The  total  cross  section  Q  is   given  by 
(3M  Q    =    rQ+rQ+     * 


rr  ^s^  IT  ^t 

where  Q  and  Q.  are  the  singlet  and  triplet  cross  section  respectively. 


(35) 


At  zero  incident  electron  energy  we  have 


Una2 


The  q\iantities  a  and  a.  are  called  the  scattering  lengths  and  are 
s     t 

defined  by  the  relation 

(36)  lim   (k  cot  6)  «=  -  i  . 

k  ->  C  * 

Using  the  data  of  Tables  I  and  II  one  can  easily  calculate  the  total 
cross  section  at  zero  incident  electron  energy.  The  results  sire  given  in 
Table  III: 

Table  III 


l«s 

^. 

Q 

Without  Polarization 

66.51^  na^ 

0 

16.61  na^ 

0 

83.15  Tia^ 
o 

With    Polarization 

60.00  Ra^ 

0 

16.61  na^ 

0 

76.61  na^ 
o 

From  Table  III  it  may  be  seen  that  the  triplet  scattering  accounts  for 
only  about  20%  of  the  total  cross  section  at  zero  energy.  Thtis,  the  singlet 
scattering  is  the  dominant  factor  for  the  process.  We  also  note  that  the  in- 
clusion of  the  r, „-  dependence  in  the  trial  wave  function  lowers  the  total 
cross  section  by  about  10^. 
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3»  Dependence  of  the  scattering  on  the  energy  of  the  incident  electron 

A,  Effective  range  theory.  It  can  be  shown  for  two-body  scattering  at  low 

T2]  hi 

energy  that  the  phase  shift  is  given  by  the  relation'- -J*  ^  J 

(37)  X     =     kcot5     =     -i  +  ir     k^ 

a       ^     0 

where  a  is  the  scattering  length  and  r  is  called  the  effective  range.  The 
above  resxilt  is  often  called  the  shape-independent  formula. 

We  shall  now  show  that  (37)  is  true  for  the  case  of  three-body 
scattering. 

We  consider  the  wave  equation  for  a  particle  of  wave  number  k  incident 
on  a  two-body  system  and  the  wave  equation  for  a  zero-energy  particle  incident 
on  the  same  tvjo-body  system.  They  can  be  written  in  the  form 

(38)  [7^+  V2  -  ^  *  ("^  *  ^^M^  "  °  * 

(39)  [7^+  V^  -  W  -  b]  T^ 


0 


where 


W  -  i^V   and   B  =-  ^  |Eg(  i 
n  h 


where  E  is  the  binding  energy  of  the  two-body  system, 
B 

We  also  consider  the  asymptotic  forms  of  the  wave  equations  (38)  and 
(39): 
(ho)  fv^+  7^+  k^-  b|  I   =  0  , 

(Ul)        jv^+  v^-  bJ  1^  =  0  . 

The  solutions  of  equations  (hO)  and  (I4I)  can  be  written  as 
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(1.2) 


I 


1 


'sin(kr,+  5)  sin(kr-+  5) 

T-    fn^'^2^  i      r.,  sin  6        '•n'-'l' 


r,    sin 


Tp  sin  6 


(h3) 


i- 


/Bn 


^  "   -  X'  ^„(^2'     1    FJ  (l  -  t)   ^n  (■■2) 


where  (^  (r)  is  the  normalized  wave  function  of  the  two-particle  system.  The 


n 


wave  function  T  and  Y  must  be  well-behaved  at  the  origin, viz: 


lim   f  <  00   and    lim 
r  -»  0  r  -»  0 


l/vi] 


(iiU) 


If  we  now  multiply  (38)  by  Y  and  (39)  by  Y  and  subtract  we  obtain 


Similarly  from  equations  (U2)  and  (I}3)  we  get 


(l45) 


,2.  „2 


2,  „2, 


|^(V^*  Slpl   -  f(V^+  S!pJ^-^  ^     II^=     0 


we  now  subtract  {kU)   from  (U5)   and  integrate  over  r,   and  r2  and  apply  Green's 
theorem.     We  obtain  after  rearranging  terms 


(U6) 


-  I    <^'^l[|    ^22[W  -  ^2^o]   *    j    ^4foV  -  ^2^o]] 


+    k 
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where  S  is  a  large  sphere  and  S  is  a  small  sphere  near  the  origin  (see 
Figure  2),  both  with  centers  at  the  origin. 


Figure  2 
As  we  let  r  become  very  large,  I  ->  J  and  ¥  -^  ^,  so  that  the  integrals  over 
the  large  sphere  <;ancel.  Furthennore  if  we  let  the  radius  of  the  small  sphere 
approach  zero  the  integrals  of  ¥  and  ¥  also  approach  zero  because  of  the  con- 
ditions imposed  on  ¥  and  VY  at  the  origin.  Equation  (U6)  thus  reduces  to 


(li7) 
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If  we  substitute  equations  (1(2)  and  (Ii3)  into  ihl)   and  carry  out 
the  integrations  v,-e  obtain 

(h8)         -  k  cot  5  -  i  +  k^  [  d  Z^d  r^   \\  \^-   Y  Y  1  =  0  . 

We  can  expand  the  wave  functions  ¥  and  ^  as  functions  of  the  wave 


nixraber  k: 

o      1 


{U9)  I  =  Y  +  Y^k  +  Y„k^  +  ...   , 


(50)  I  "  lo*  ^1^  "^  ^2*^^  "^  •••   » 

2 

so  that  (It8)  can  be  written  to  order  k  as 

(51)  k  cot  6  -  -  i  ""  ^^  I  ^^i'^  ^2  r^"  ^ol  • 
If  we  define  the  effect! -^/e  I'ange  r  as 

then  (5l)  becomes 

(53)         k  cot  5  =  -  I  +  ^  ^Q^^     > 

which  is  the  effective  range  forriiula. 

If  we  are  dealing  with  a  bound  state  of  the  three-particle  system 

Me   can  use  the  same  arguments  as  were  used  in  obtaining  (53),  except  that  the 

asymptotic  wave  function  J  must  now  vanish  for  very  large  r.  The  result 

of  carrying  out  this  argument  is  to  obtain  a  relation  between  the  attachment 

2 
energy  -k-    for  the  third  particle,  the  scattering  length, and  the  effective 

rangej  this  relation  is 

(5U)  e  -  I  .  I  r^  e^  . 
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For  S-wave  scattering  the  cross  section  can  be  written  as 

(55)  Q  =  ^  sin^  6  . 

k 

Kor  low  energies  we  can  use  (53),  so  that  (55)  becomes 

(56)  Q  =   g   ^"  ^   -        ^" 


k  +  k  cot  6     ,  2 
k  + 


[-M^^^I 


We  can  thus  see  that  the  total  cross  section  (singlet  plus  triplet) 

is  determined  by  the  four  parameters  a  ,  r  >  &x  and  r  .  • 

We  ha-ve  already  calculated  a   and  a  .  .  Equation  (56)  is  valid  so 

OS      ot 
2 
long  as  (r  k)  is  negligible  compared  to  (r  k).  We  shall  see  that  the  limit 

of  applicability  is  determined  by  r  •  For  our  problem  we  can  use  (56)  for 

OS 

values  of  k  up  to  O.Oij  or  about  0.C2  ev.     This   is  of  the   order  of  thermal  energy 

We  now  seek  a  method  of  calc\ilating  r       and  r  ,  •     Since  a  bound  state 

^  OS      ot 

(singlet  state)  exists  for  the  k"  ion  we  can  use   (5Ii)  to  find  r   if  we  know 

OS 

the  attachment  energy  of  the  second  electron.  This  quantity  has  been  calculated 

by  Chandr£sd<faar  (result;  0.7019  volts)  and  measured  by  Bransccmb  I- -'  (result 

0.756  volts),  ^or  the  triplet  scattering  no  corresponding  bound  state  exists, 

so  that  we  shall  have  to  use  other  means  to  determine  r  . . 

ot 

Returning  now  to  the  calculation  of  r  ,  we  shall  vise  the  experimental 
value  of  the  wave  number  of  the  second  electron,  which  is  £  =  0.236.  Substitut- 
ing this  value  of  e  and  the  values  of  a  into  (5U)  we  find  the  results  shown  in 

s 

Table  IV: 

Table  IV 

r   (without  polarization)     r  (with  polarization) 

r^„       it.GU  a^  3.85  a 

OS  o  o 
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B.  Independent  Calculation  of  r  .  The  quantity  r  .  must  be  evaluated  by  a 

'      0  ox 

method  independent  of  the  effective  range  theory  since,  as  mentioned  above,  no 
triplet  bound  state  of  H"  exists.  If  we  choose  a  trial  function  which  depends 
on  the  energy  of  the  incident  electron,  then  (25)  and  (28)  will  give  X  =  k  cot  5 
as  a  function  of  the  energy  of  the  incident  electron.  In  (28)  the  quantities 
a,  b,  ...  and  A  and  B  will  also  depend  upon  the  energy  of  the  incident  electron. 

For  electrons  whose  energies  are  close  to  zero,  we  can  expand  each 
term  in  (28)  in  powers  of  k.  According  to  (53),  we  obtain 


(57) 


X  «  Pq+  P-j_k^+  O(k^)  . 


Differentiating  X  with  respect  to  k  ,  we  have 


(58) 


(59) 


ax 
ak^ 


Pj+  o(k)  . 


In  the  limit  k  -^  0,  we  have,  using  (37) > 


k  -^  0  Vaky 


V 


Pi  =  ?  ^c 


We  now  choose  an  energy-dependent  trial  function  having  the  correct 
asymptotic  form,  namely 


(60) 


h  '    ^t 


6     e     +e     8 


^ 


-rr. 


V^n  L     1 


-rr,  2 

(1-e    ) 

, sin  kr,  +  -^ i—  cos  krT 

k  r,         1        r,  1 


B(l-e   ^) 


e    +  r^ 
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Then  equation  (25)  can  be  written  as 

(61)  X     •=     -   [aA^+  bB^+  cAB  +  d+eA  +  (f-DBj      . 


2 
Let  us  now  find  3X/3k  .     We  have 


/^„x  ax  ^a     .2     db     ^2     dc     .r,  ^  3d         de     .   ^  Bf    „ 

(62)  -  -=5     =     — «•  A  +  — X  B  +  — J  AB  +  — K  +  — J'  A  +  — k  B 

dk  ak  ak  ak  ak       ak  ak 


.  2aA  ^  .  2bB  ^  .  c  f  A  ^  .  B  ^)     .  e  ^  .   (f-1)  ^ 

ak  ak         Vs.   ak  aky  ak^  ak 

Let 

(63)  a     =     a  +  a-k  +  a^k  +  ,..,         b     =■     b  +  b-k  +  'bjr+  ,..,     etc. 


so  that  as  k  — *0,   equation  (62)   becomes 

(6h)  ^lim     f^    -     a^A^.  b^B^.  c^^A^B^-  d^.  e^A^.  f^B 


ak 


r2a  A  +  c  B  +  e  If^  )  +    r2b  B  +  c  A  +(f-l) 
00       00       oK^,  2/       loo       ooo 


'vak, 

o 


From  the  conditions  (29)  and  (30)  we  see  that  the  last  two  terms  in 
(6U)  vanish,  so  that  we  obtain 

^^^^      J-"^   n    fe)   =   I  ^O  •=   -  Fi^O^  ^1^0*  ^1^^*  ^'    ^1^0^  Vol   • 

k  — >  0  \ak/  •-  -J 

In  the  numerical  calculation  of  lim   ( — ?r)  we  shall  use  the  value 

k  -»  0  \^ak  / 
of  Y  which  gave  us  an  extreinum  for  X  .  This  assumption  may  be  justified  by  the 

following  argument.  We  may  write 

(66)  X  =  f(A,B,Y)  +  g(A,B,Y)k^  , 

where  A,  B,  and  y  ^.re  functions  of  k. 
Then 


(67) 


-  22a  - 
dk  dk  dk^ 


dA      ^   dr  db      _^   dl    dv     ^      / .    _      \    ^  dg     ,  i 

-T  "  3B  T?     a^  f?     ^^^'^''^^     t^  ^ 

dk  dk  '    3k  dk 


Letting  k  ->  0  gives 

where  the  subscript  o  indicates  the  values  A  =  A  ,  B  =  B  and  Y  =  Y  »  for  k  =  0. 

From  the  fact  that  X  is  an  extremum  for  A  =  A  ,  B  =  3  and  y=  Y  it 

o  o'     o       'o 

follows  that 

Wo '  (-1 '  (i  =  ° 

so  that   (68)   reduces  to 

(70)  lijn      {  ^  ]    =     g(A_,B_,Y J     =     P.    -     i  r 


.^0©    °     '''o'^^"o^     '     "1-    1-0 


All  of  the  integrals  involved  can  be  evaluated  by  elementary  methods 
(see  Appendix  III). 

The  results  are  shown  in  Table  V: 

Table  V 

Singlet         Triplet 
r         3.11  a         0.808  a 

0  0  0 
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The  above  value  of  r   differs  from  that  calculated  usinc  the  effective 

OS  ° 

range  theory  by  about  30/6.  The  present  value  for  r   is  probably  not  as  ac- 

OS 

curate  as  the  one  obtained  from  the  effective  range  theory  siix;e  the  trial  func- 
tion vra,s  not  as  accurate. 

In  Section  U   we  shall  discuss  the  sensitivity  of  the  cross  section  to 

the  parameters  a  ,  r  ,  a.  and  r  .  .  V.^o  sfjall  show  there  that  the  low-energy 
S   OS   t      ot 

cross  section  is  not  appreciably  influenced  by  errors  in  r  . .  Hence  the  above 
method  for  calculating  r  .  need  not  be  discarded  because  of  its  inaccuracy. 
For  the  singlet  cross  section  we  shall  use  the  value  of  r   and  a 

^  OS       s 

calculated  using  the  polarization  approximation.  For  the  triplet  cross  section  we 
shall  use  the  value  of  a  calcxilated  for  the  non-polarization  approximation, 
and  the  value  of  r  .  found  from  the  energy-dependent  trial  f\inction. 
The  results  of  this  calculation  are  given  in  Table  VI: 


Table  VI 

k 

^s/-c 

^t/-o 

k     s'      0 

^t/-o 

Q/na^ 

0.000 

2U0.0 

22. lU 

60.00 

16.61 

76.61 

0.005 

2lj0.0 

22.U 

60.00 

16.61 

76.61 

0.010 

239.2 

22.12 

59.81 

16.59 

76.I4O 

0.015 

238.2 

22.11 

59.56 

16.58 

76.1)4 

0.020 

237.2 

22.10 

59.31 

16.58 

75.80 

0.025 

235oh 

22.08 

58.56 

16.56 

75.ij2 

0.030 

233.5 

22.05 

58.38 

16. 5U 

7U.92 

0.035 

231.1 

22.02 

57.77 

16.52 

71.29 

o.oho 

228.6 

21.96 

57. U 

I6.I47 

73.61 

These  results  are  plotted  on  Graph  I. 
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l4.  Sensitivity  of  the  method 

For  three-body  collisions  there  is  no  rigorous  way  of  estimating 

upper  and  lower  bounds  for  the  phase  shifts,  so  that  we  have  no  means  of 

detennining  the  accuracy  of  our  calcu!lation. 

We  shall,  however,  determine  the  effect  which  errors  in  a  and  r 
■»        '  o 

have  upon  the  cross  section  for  the  scattering  process. 

The  total  cross  section  as  given  by  equation  (31)  is 

(3U)        '  ^    =  i%'  ^%    ' 

Suppose  there  is  an  error  ^  Q  and  ^  Q.  in  tiie  values  of  Q  and 
Q.  respectively;  then  Q  will  be  changed  by  an  amount  ^  Q.  The  relative  error 
^  is  then  given  by  the  relation 

(71)  ^  "   ^ ^—   . 

It  can  now  be  shown  that  an  error  -— -  in  a  and  in  r^  leads 

a  r       o 

AQ  ° 

to  an  error  -^    in  Q  for  either  Q  or  Q.  j  this  is  given  by 

(72)  ^    . SJ-Sf — -J    -    1    . 


k^ 


Ar 


(73) 


For  k  =  0  equation  (72)   becomes 
AQ     .      o  Aa  /^a\2 


We  make  the  supposition  that  the  trial  fiinction  for  the  singlet 
scattering  is  a  good  one  so  that  the  value  cfa  and  r      which  we  have  calculated 

S  OS 

are  quite  acciirateo 
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For  the  triplet  scattering  such  accuracy  is  very  unlikely,  as  pointed 

out  previously.  We  shall  therefore  assume  that  any  error  in  the  cross  section 

is  due  to  errors  in  a.  and  r  . • 

t     ot 

Before  proceeding  to  examine  the  effects  that  errors  in  a.  and  r  . 
produce  in  the  cross  section  let  us  compare  our  results  for  zero-energy  scat- 
tering with  those  obtained  by  extrapolating  the  calculations  of  I^ssey  and 

ri9i 

Moiseiwitsch"-  -■  to  zero  energy.  We  shall  use  the  values  for  what  the  above 

authors  call  the  'exchange-polarization  approximation.' 

2 

We  carry  out  the  extrapolation  by  plotting  k  cot  5  against  k  (see 

Graph  II) •  For  both  the  singlet  and  the  triplet  case  the  points  fall  almost 
on  a  straight  line.  The  results  obtained  from  this  graph  are  shown,  along  with 
our  results,  in  Table  VII: 

Table  VII 

k  cot  5  r 

o 

M  and  M         Our  Value         M  and  M  Our  Value 

Singlet  -O.lhS  -0.1291  2,67  a  3.85  a 

°  0  0 

Triplet     -0.1i3     -0.1^251     1.56  a^     0.308  a^ 

Our  values  for  r   and  r  .  differ  quite  markedly  from  those  of  Massey 

and  Moiseiwitsch.  The  values  of  r  obtained  by  extrapolation  and  by  calculat- 

2 

ing  lim   (dX/3k  )  are  probably  not  as  good  as  those  obtained  from  the  effective 

k  ->  0 
range  formula,  since  the  former  method  does  not  insure  that  a  good  value  of 

the  intercept  will  yield  a  good  value  for  the  slope  of  the  line,  while  the  latter 

method  is  founded  on  strong  theoretical  grounds. 
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Moreover,  it  can  be  seen  from  (72)  that  the  cross  section  is  in- 

2 
sensitive  to  r  for  small  values  of  k  ,  so  that  even  relatively  large  errors 

in  r  will  not  have  an  appreciable  effect  on  the  cross  section. 

Our  value  for  a  differs  from  the  value  obtained  by  extrapolation 
s 

from  the  results  of  Massey  and  Moiseiwitsch  by  about  iVo  and  should  be  re- 
garded an  an  jjnprovement  due  to  the  use  of  our  more  acciirate  trial  function. 

The  values  for  a  are  about  the  same  for  both  calculations. 
t 

The  zero-energy  cross  section  obtained  by  extrapolation  from  the 

2 

Massey  and  Moiseiwitsch  results  is  6h.6  na  ,  as  compared  to  our  value  of 

76.61  na  ,  which  is  about  20^o  higher.  No  experiments  heve  been  done  in  the 
energy  range  where  our  calculations  are  valid,  so  that  there  is  no  way  of  com- 
paring our  calculation  of  the  cross  section  with  experiment.  Recent  experi- 
ments L  -I  in  the  region  between  two  and  ten  volts  give  values  of  the  cross  section 
higher  than  those  calculated  by  Massey  and  Moiseiwitsch  for  these  energies. 
It  thus  appears  likely  that  our  larger  zero-energy  cross  section  is  an  improve- 
ment over  the  extrapolated  calculation  of  Massey  and  Moiseiwitsch. 

Let  us  now  examine  the  effect  of  errors  in  a.  and  r  .  on  the  calculated 
value  of  the  cross  section.  As  mentioned  above,-we  take  the  point  of  view  that 

a  and  r   are  accurately  known  and  that  a.  and  r  .  may  be  in  error, 
s      OS  t      ot   *' 

Let  us  assume  that  there  are  large  errors  in  both  a.  and  r  . .  We 

shall,  for  the  sake  of  argument,  take  -—^  «=  0.50  and  — =  1.00.  If  we 

^t  ot 

now  substitute  these  values  in  (72)  and  apply  (71)  we  obtain  the  following  values 

for  the  relative  error  in  the  cross  section  produced  by  the  errors  in  a^  and  r^^J 
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Table  VIII 


-»—  (relative)         -w—  (percentage) 


O.OCO  f§lg  27.1 

0.005  fg^  27.0 

0.010  fg^  27.1 

0.015  '^  27.2 

0.020  ^^  27.3 

0.025  1^  27.1^ 

0.030  f^  27.5 

0.035  f^  27.6 

O.OhO  §lg  27.7 


Prom  "feble  VIII  we  see  that  the  asstimed  errors  in  a.  and  r  .  produce 
errors  in  the  cross  section  which  are  less  than  30 /b.  These  results  are  shown 
on  Graph  III.  If  the  errors  in  a.  and  r  .  are  no  larger  than  we  have  supposed, 
then  the  true  cross  section  will  lie  between  curves  A  and  B  on  Graph  III. 

Until  more  rigorous  methods  are  obtained  for  estimating  the  errors 
in  this  type  of  calciilation  the  above  estimate  will  have  to  suffice. 
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5»  Neutral  impurity  scattering  in  semiconductors 

Measurement  of  the  electrical  properties  of  germaniiira  by  Lark- 
Horowitz  and  Johnson'-  -^  and  silicon  by  Pearson  and  Bardeen^  J  show  that  the 
temperature  dependence  of  the  resistivity  of  these  materials  can  not  be  ex- 
plained by  lattice  scattering  alone,  especially  at  low  temperature.  In  order 
to  account  for  the  resistivity,  the  effect  of  scattering  by  ionized  and  neutral 
impurities  must  be  considered. 

Conwell  and  Weisskopf  •-  -^  have  calculated  the  effect  of  ionized  im- 
purity scattering  using  the  Bom  approximation,  and  Erginsoy^  -^ ,   using  the 
results  of  Massey  and  MoiseiwitschL  -^,  has  calculated  the  effect  of  neutral 
impurity  scattering. 

The  scattering  of  holes,  or  electrons,  by  neutral  impurities  in  serai- 
conductors  is  similar  to  the  scattering  of  electrons  by  hydrogen  atoms.  In  a 
seiJiiconductcr  the  radius  of  the  first  Bohr  orbit  of  the  electron,  or  the  hole, 

is  increased  by  a  factor  D,  the  effective  dielectric  constant  of  the  material, 

2 

This  would  lead  to  an  increase  in  the  cross  section  by  a  factor  of  D  . 

We  shall  now  obtain  an  expression  for  the  resistivity  as  a  function 
of  the  cross  section. 

The  resistivity  can  be  written  as 

C7U)  P„-  -^  , 

where  m  is  the  effective  mass  of  the  conduction  electron  or  hole,  e  is  the 
charge  on  the  conductacn  electron,  n  is  the  concentration  of  conduction  electrons 
or  holes.  T  is  called  the  relaxation  or  collision  time  and  is  related  to  the 
cross  section  by 

(75)  T  =  -V  ' 
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where  N  is  the  concentration  of  neutral  impurities  in  the  semiconductor. 
Combining  equations  (7U)  and  (7$)  gives  as  our  expression  for  p 

kffi  \ 


(77) 


n 


e.g. 3 


N 


=  3.65  X  10'^  kQ  —  ohm-era 

n 


vrhere  k  and  Q  are  in  atomic  units,  and  D  =  13» 

V^e  have  calculated  p  for  silicon-boron,  sample  2  and,  silicon- 

^  r2ii 

phosphorous,  sample  B,  which  were  investigated  by  Pearson  and  Bardeen'-  -' . 
The  results  are  shown  and  compared  with  those  of  Erginsoy  in  Tables  X  and  XI. 

Table  X 
p  (Silicon-Boron,  Sample  2) 

Erginsoy 


lA  (°K)-i 

p   (ohm-cm) 

lA  (°K)-^ 

p  (ohm-cm) 

p, (ohm-cm) 

0.01018 

0.0521 

0.012 

0.0866 

0.0526 

0.007072 

0.0177 

0.010 

O.GI4O2 

0.0390 

0.005198 

0.00725 

0.006 

0.0086 

0.0191 

0.003978 

0.00U6h 

0.005 

Table  XI 

0.00h6 

0.0153 

p  (Silicon' 

-Phosphorous, 

Sample  B) 

Erginsoy 

lA  (°K)-i 

p  (ohm-cm) 

lA  (°K)-i 

p  (ohm-cm) 
n 

p.  ( ohm-cm) 

0.01018 

0.0103 

0.012 

0.0077 

0.0283 

0.007072 

0.00862 

0.010 

0.0077 

0.0250 

0.005198 

0.00U81 

0.006 

0.0039 

0.0155 

0.003978 

0.002142 

0.005 

0.0029 

0.0129 

.p.  is  the  resistivity  due  to  scattering  by  ionized  impurities. 
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The  results  are  also  shown  on  Graphs  IV  and  V. 

Our  results  for  the  resistivity  are  larger  than  those  found  by 
trginsoy,  which  is  to  be  expected  since  our  values  for  the  cross  section  are 
larger  than  those  of  Massey  and  Moiseiwitsch. 

These  results  indicate  that  the  effect  of  neutral  impurity  scatter- 
ing in  semiconductors  is  of  more  Importance  than  had  previously  been  thought. 
This  is  especially  true  at  very  low  temperatures  where  the  effect  of  lattice 
scattering  is  negligible. 
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Appendix  I.  Calculation  of  variational  parameters  using  T.  =  ^  ■*"  /f 

t 

We  consider  the  integral 


/  ^t^  ^t  ^  ^1 


d'i',   , 


where  Y.  is  given  by  equation  (2?)  and  L  is 

(I.l)  L  =  v2^  v2  +  2(i.  +  A-  -  ji-  -  |)  . 

±       d  ^1   ^2    12  '^ 

The  terms  a  ,  b   ,  etc.  in  (b)  are  then  obtained  by  evaluating 


the  follovdng  integrals: 


f  r"*^i  "^^p    "P^i  •"^2]     r~^i  "p^2    "^^1  "^2! 

1^  -  e       -^  e       "^  +  e       ^  e       "^     X    L  e       "^  e       '^  +  e       "^  e       "^  d  r^d  t^  j 

'0=    ^/    [(1-"^V2  i  (l-e"'"^)e"^] 

r       -yr,     -r  -Yr„     -r"! 

X   L   l(l-e       ^)e     ^  +  (1-e       ^)e     "M  d  r^ 


dr,  , 


1 
^5"  n    [_' 


j     e  e  +  e  e 

L     (1-e       ■^)e     "^  +  (1-e       ^)e     "^     d  r j_ 

I  [(l-e"^l)e"''2  +  (l-e"^2)g-^lj 

r-or,     -prp  -fr       xir  1 

Le       -^e       '^  +  e       "^e       ^d^^dr 


dr. 


i7 


X  L 


(1-e       h  e     ^     ^     (1-e       ^)  e     ^ 


-YTt   2  -r,  -^<,  2  -r, " 

(1-e      •^)  e     ^    ^     (1-e      ^)  e     "^ 


dZ^  dt^  J 
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1 


-YTt   2  -r 


j    [_e  e  .  e  e        J 


-YTp  2  -r - 
)  e     ^     ^     (1-e       ^)  e     ^ 


^V^2 


(1-6       ^)  e     ^       ^     (1-e       ^)  e     "^ 


-ar,     -pr„        -pr. 


[;-""i 


e  +  e 


'^T^^i 


df. 


— ? 


[/  \l-^' 


-r. 


-rr2     -r 


(1-e      ^ 
.  "1 


)e     ^  +  (1-e       ')e 
-YT^   2  -rg 


'] 


-YT    2  -r,"! 
)   e     '     ^     (1-e       ^)  e     ^ 


dr^  d^^ 


(1-e      h   e     ^     ,    (1-e       ^)  e    ^ 


(1-e       -^)e     ^  +  (1-e       ^)e 


i]ar^ 


df. 


In  order  to  carry  out  the  above  integrations  we  raiist  find  the  effect 
of  L  on  the  functions  upon  which  it  operates.  The  res\ilts  of  operating  with  L 
are  listed  below: 


>,  p2_  ,  ,  2(1^  ,  2(i:£)         2  1     -1 


Lie     ^^e"^"^] 


,2,  p2.  ^  ^  2(i:i)  ,  2(i:a)   .  ^n     ^^1 
^1  ^^2  ^12  J 


-pr,     -or 


2 


-  33  - 


[(l^-^^)3-^^]        .      ^ 


-YTi     -^2 


2     ^^1       ^2^/2  2   w^       ^Iv     ^2 

■- re  e         +  (-r  -  r— )(l-e         )e       J 

^1  U       ^12 


.[( 


-YT      -r^ 


e 


2     -^2     """l       ,2  2  w,     "^2,   "^1 


LHl-e       ')e 


fj     „     ^r^     ^^  *- y'  e       ^e     "M^-^)(l-e       ')e     "^j 


■2       ^2 


(1-e       •^)  e     ^ 


hY  e         ^e     ^       2rS       ^e     ^  ^  .i_    J_.    (l-e       ^)S     ^   . 
^1  ^1  ^1     ^12  ^1 


r        -YTp  2  -r  ■ 


Uy  e        ^e    -^  _  2y  e       ^e     "^  ^  (2..  -L)   d-e         )  e    "^  ^ 
^2  ^2  ^2     ^12  ''2 


Using  the  above  relations  the  integrals  can  be  evaluated  by  elementary- 
methods.     The  results  of  these  calculations  give  for  a  ,  b  ,  etc.: 

°^  o       o 


a     =     32n 
o 


2.  „2 


^^^(^^-4f^-u..,.^l 


b     •    U 


1         1  Y^-l4  Y*-2 

,^  "  ^  "  (y*2)^  "  8(Y+l-)l 


Jni 


(1+y) 


^E'H^"^^H^] 


+  u 


p- 


Y^-3  2y+3 


(y+2)- 


,     3(y^1) 

(y+1)^(y+2)^       8(y^1)^ 


0" 


12  8n 

>5 


1       f       ay       ^   Q+2(l^^^)     _     ajxtiO^f) 
(l^  Ua+Y)^         (a+p+1)^  (a-^^+Y+1)^ 

+        1         /py         +  p-^2(l+a)     _     pn-+2(l+a/ 


-  3U- 


d     = 
o 


^2  log(l+  J)   -  3  log(l+Y)   +  2  logd  +  ^  )   -  ^  log(l+2Y) 
.1         2      .        3  2  1 A 


^iiii 


il 


2nr 


(h-2y)^     (in)' 


1  2        ^  1_^ 

L     On?     (i+2y)' 


+  2 


1         2  1 "I 

?-^^?n:TYTj 


ii^ril 


1     A.  +  1 


(1+2y)' 


^Tl^-^^^TiWJ 


6Un 
/2 


2y^ 


(1-t^)^    L(^+2y) 


■7  -     ^      .1   *  -^  [20130^ 
)2       (an)^J        (1-^)^    L(a-^+l) 


_     2[2(an)^3(l^)]     ^     2(a.2Y)^3(l^) 


(a-»f+Y+l) 


(a-.f+2Y^l)' 


ix 


2     I- 


.(1+a)- 


(P+2y)         (P+y) 


2[2(p+Y)+3(l+a)]  2(B 


(a+p+Y+1)' 


]) 

■»-2Y)+3(l+a)"1 
.B+2y+1)^        J 


3(l+a) 

"7F 


f    -    h 


19      1 


+  ii 


V3       _     2y+3       +      Y-H 


(y+2)^     ii(Y+l)  (3Y+2)' 


.(1+t)- 


[•• 


?  * 

(1+y)  (1+2y) 


W-- 


(in)^(2+Y)^ 


S^r;:br 


h(l+Y)^l+2Y)^  (2+y)^       2(l+r)^  (1+2y)^(2+3y)^ 


Natviral  log. 
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Appendix  II.  CalciLlation  of  variational  parameters  using  T,=  (l+cr, p)  J  +  X, 
We  consider  here  the  evaluation  of  the  integral 


/  ^"t'*'^! 


dr^ 


where  Y.  is  given  by 


(33)         Y^  =  (1  +  c^r^2^  1  +  Z  . 
We  can  then  write 

I     =     P+Q+R  +  S, 


where 


p  =  j  (l+Z)  Lil*X)  dr^  dr^ 

Q  -  f  (|+%)   L(c^r^2  1)    ^"^1  d'i'2 

R  =  J  C3_r^2  1  ^^I  +?^)  d  r^^  d  r^ 

^  "  J  °1^12  5  ^^*^l'"l2  ^  d  r^^  d  ^2 


P  has  already  been  evaluated  in  Appendix  I, 

It  can  then  be  easily  shown  by  using  Green's  theorem  that  Q  «  R« 
Since  we  have  already  evaluated  the  quantity  L(^  +%)   it  id.ll  be  simpler  to 
\xse  the  integral  R. 

The  integral  S  can  be  simplified  by  putting  in  the  explicit  value  of 
L.  We  can  then  write 

S  =  S^+  Sg  , 


where 


Si  ■  =1  I  ''12  1  <'l*  '2"''l2  i  '  """l  *'^2  ' 
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Sp  can  not  be  simplified  any  further. 

Fron  the  symmetry  of  ^  it  is  clear  that  we  have  ^t  J  "  ^  $•  3- 
can  then  be  written  as 

h"     ^4    j    ^12i^l  ^^12!^  ^'^1^^2  • 
It  is  then  not  difficult  to  show  that  S,  takes  the  form 

v^ere  7.  is  the  usual  gradient  operator  and  vrhere  we  have  used  the  fact  the 

2 

(V-  r,  p)  =  1«  All  of  the  above  integrals  can  again  be  evaluated  by  elemen- 
tary methods. 

The  only  terns  which  will  be  affected  by  the  r^^-term  are  those  which 
depend  on  a  and  pj  these  are  a  ,  c  and  e  • 

Below  we  list  the  additional  terms  a  ,  c  and  e   which  must  be  added 

o'     o  0 

to  those  calcxilated  in  Appendix  I  for  a  ,   c     and  e   : 

^^  o*     o  0 


t  2 

a     -    32n  c, 
o  1 


2     .2     ,    I- 


g  +  p  -  1 
I6(a^)^ 


UO         .         20         .         9 


.   3    .      ho 


p(a+pr       p    {a+pr       P    (a+p)       ^       a(a+p) 


a   (a+p)  a  ia+p) 


■i] 


+     (l-<^)  r     16        ^      10  6        ^  3    ^ 16_  ^        6        ^        2     "1 

8(a+p)^  [_p(a+p)^       p^(a+p)^     p^(a+p)       p^      a(a4p)3     a^(a+p)^  a^(a+p)J 

^     (1-P)  r    16        ^      10         ,         6        ^   3    ^ 16_  ^        6        ^        2     "I 

8U^^  |_IU+P^       a^Ca+p)^      a^(a+p)       ?       p(a4p)^     p^(a-^)^  p^(a+p)J 


_1     ^   ,    f  3$(a^^  p^-  1)     ,     bh  -  30(a-f ) 
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8nc, 


1       2' 


r 


8(a+p)- 


UP       ^       2C       ^       9       ^  3 
p(a+p)^  p^(c+p)2  B^(a+p)  ^ 


1         r  -^  1  3       p.    ^   iH  218    _      96      _      1      ^      16 


8(a+p; 


38l4ap     _     3 
(o+pT       hap- 


Fi  ^  i~l     3   n 

[a         pj        Ua-^p    [a 


^    +    1 

a         p 


32lTC^ 


8y 


(a+p+Y+1)' 


16    __      ^  10  ^  6 

(l+p)(a+p+r+a)-^       (l+plF(a+p+Y+l)^       (l+p)^(  a+p+Y+1) 


+       3       +  16  ^  6  ^  2  "1 

(1+^)^        (a+Y)(a+fi+Y+l)  (a+Y? (a+p+Y+1)  ( a+Y)^+p +Y+1 )    I 


sy 


2         r- 


( a+p+Y+1) 


5 


Uq ^  20  _^  9 

(l+p){ a+p+Y+1)^        (l+f)^(a+p+Y+l)^        (l+p)^( a+p+Y+1) 


+       3  ho  ,  20  9  ^3 

(l+p)^   (a+Y)(a+p+T*l)      (a+Y)    (a+p+Y+1)        (a+Y)^(a+p+Y^■l)    (o+y)^ 

+  8         r  lb  ^  10  ^  6  ^       3 

(a+p+1)^    |_(l+p)(a+p+l)^       (l+p)^(a+p+l)^        (l+p)3(a+p+l)        (l+p)^ 
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Appendix  III,  Calculations  required  for  r   and  r  .  '/dth  energy-dependent 

OS      ox  o^    i 

trial  functions 
"We  wish  to  calculate  the  quantities  a,,  b, ,  etc.  in  equation  (63) » 

In  order  to  accomplish  this  we  must  evaluate 
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where  f,  is  the  energy-dependent  trial  function  given  in  (60)  and  L  is 

L  =  V^*   V^(-^  +  i_  »  -i_  .  1  +  k^)  . 
1   2  r^   r^   r^2 

To  carry  out  the  calculation  we  must  find  the  result  of  operating 
on  the  asymptotic  part  of  (63)  with  L.  These  results  are  listed  below: 
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The  results  for  r,  and  r_  interchanged  are  exactly  the  same  as  those 
above  with  r,  and  r^  interchanged* 

With  these  results  we  can  now  carry  out  the  calculation  of 
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We  now  list  the  values  of  a^ ,  b^ ,  etc. 
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